Controlling domain patterns far from equilibrium 
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A high degree of control over the structure and dynamics of domain patterns in nonequilibrium sys- 
tems can be achieved by applying nonuniform external fields near parity breaking front bifurcations. 
An external field with a linear spatial profile stabilizes a propagating front at a fixed position or in- 
duces oscillations with frequency that scales like the square root of the field gradient. Nonmonotonic 
profiles produce a variety of patterns with controllable wavelengths, domain sizes, and frequencies 
and phases of oscillations. 



PACS numbers: 05.45. +r, 82.20Mj 

Technological applications of pattern forming systems 
are largely unexplored. The few applications that have 
been pursued, however, have had enormous technological 
impacts. Magnetic domain patterns in memory devices 
provide an excellent example [1] . Intensive research effort 
has been devoted recently to dissipative systems held far 
from thermal equilibrium [2]. Unlike magnetic materials, 
such systems are nongradient in general and their asymp- 
totic behaviors need not be stationary; a variety of dy- 
namical behaviors can be realized, including planar and 
circular traveling waves, rotating spiral waves, breath- 
ing structures and spatiotemporal chaos. This wealth 
of behaviors opens up new opportunities for potential 
technological applications. Their realizations, however, 
depend on the ability to control spatiotemporal patterns 
by weak external forces. Most studies in this direction 
have focused on drifting localized structures [3]. 

In this paper we present a novel way to control do- 
main patterns far from equilibrium. We consider dissipa- 
tive systems exhibiting parity breaking front bifurcations 
(also referred to as nonequilibrium Ising-Bloch or NIB 
transitions [4,5]), in which stationary fronts lose stability 
to pairs of counterpropagating fronts. Examples of sys- 
tems exhibiting NIB bifurcations include liquid crystals 
[6] and anisotropic ferromagnets [7] subjected to rotating 
magnetic fields, chains of coupled electrical oscillators [8], 
the catalytic CO oxidation on a platinum surface [9,10], 
the ferrocyanide-iodate-sulphite (FIS) reaction [11], and 
semiconductor etalons [12]. A prominent feature of these 
systems is that transitions between the parity broken 
states, the left and right propagating fronts, become fea- 
sible as the front bifurcation is approached. Indeed, in- 
trinsic disturbances, like front curvature and front inter- 
actions, are sufficient to induce spontaneous transitions 



and can lead to complex pattern formation phenomena 
such as breathing labyrinths, spot replication [11,13] and 
spiral turbulence [14]. It is this dynamical flexibility near 
NIB bifurcations that we wish to exploit. By forcing 
transitions between the left and right propagating fronts, 
using spatially dependent external fields, we propose to 
obtain a high degree of control on pattern behavior. 

We demonstrate this idea using a forced activator- 
inhibitor system of the form 

Ut = e~ (u - v? — v) + u xx , 

v t = u - a\v + Sv xx + h + Jv Xl (1) 

where u, the activator, and v, the inhibitor, are scalar 
real fields, and h and J are external fields. With ap- 
propriate choice of a± > 0, the system (1) has two lin- 
early stable stationary uniform solutions, an "up" state 
(u+,v+) and a "down" state (u—,V—), and front solu- 
tions connecting these states. The domain patterns to 
be considered here consist of one dimensional arrays of 
up state regions separated by down state regions. In the 
absence of the external fields, a stationary front solution, 
stable for e > e c (6), loses stability in a pitchfork bifurca- 
tion to a pair of fronts propagating in opposite directions 
at constant speed. For e/i « 1 the bifurcation point is 
given by e c = where q 2 = a± + 1/2 [14]. Activator 
inhibitor models have been used to describe some of the 
systems mentioned above [8,9,10,12]. In these systems 
h is usually a parameter that introduces an asymmetry 
between the up and down states. In the context of chem- 
ical reactions involving ionic species J may stand for an 
electric field [15]. 

Consider first the effect of a constant external field on 
the front velocity. Away from a front bifurcation the 



1 



0.5 



o 
o 

o 
> 

-+-> 
C 

o 

u 



0.0 



-0.5 

-0.05 



+ 



+ 



Ji 0.00 J a 0.05 

J 

FIG. 1. Front velocity vs. external field, J, near 
the nonequilibrium Ising- Bloch bifurcation. Parameters: 
oi = 3.0, e = 0.01, S = 2.0, h = 0. 



effect of a weak field is captured well by a linear ap- 
proximation and therefore the effect of the field is small. 
This is not the case close to a front bifurcation; the ve- 
locity - external field relation becomes multivalued (or 
hysteretic) even for a weak field as illustrated in Fig. 1 
[8,10,13,16]. This form is a generic unfolding of a pitch- 
fork bifurcation and holds for various unfolding forces 
and parameters including intrinsic disturbances like cur- 
vature [13]. We emphasize that the termination points, 
Ji and J 2l of the upper and lower branches lie close to 
J = 0. The significance is that weak external fields can 
induce transitions between the two branches, or reverse 
the direction of front propagation. 

The effect of a nonuniform external field [17] can be 
understood intuitively in the following way. Consider a 
constant h and a linear profile for J: J = —ax, where 
< a <C 1. This choice divides space into three regions 
according to the type and number of existing front so- 
lutions: (i) x > x\ = —J\/a, where J < Ji and only a 
single front corresponding to a down state invading an 
up state exists, (ii) x < x 2 = —J 2 la, where J > J 2 and 
only a single front corresponding to an up state invad- 
ing a down state exists, and (iii) x 2 < x < xi, where 
Ji < J < J 2 and both fronts coexist. This profile of J 
results in an oscillating front, with oscillations roughly 
spanning the interval x 2 < x < x\. The front propa- 
gation direction is reversed during transitions from the 
upper to the lower velocity branch at J k, J\ and from 
the lower to the upper branch at J m J 2 . Obviously, a 
variety of pattern behaviors can be induced using a non- 
monotonic J profile. For example, a single hump profile 
can induce a breathing domain. 

We now turn to a quantitative study of front dynam- 
ics and relate pattern characteristics (e.g. breathing fre- 
quency) to control parameters. We assume e « 1 and 
distinguish between an inner region, the narrow front re- 



gion, x w Xf, where u varies sharply over a distance of 
0(\fe), and outer regions, x < Xf and x > Xf, where u 
varies on the same scale as v. 

Consider first the inner region. Expressing (1) in a 
frame moving with the front, x —* r = x — Xf(t), stretch- 
ing the spatial coordinate according to z = rj\fe, and ex- 
panding u = u + v feu 1 +eu 2 + ..., v = v + v ^ev 1 +ev 2 + ..., 
we find at order unity the stationary front solution, 
mo = — tanh(z/\/2), vq = 0. At order \[e we find 



C,U\ =V\- XfUOz 

Solvability of (2) gives 



Vf = 



VV2 



C = d 2 z + l- 3u 2 . (2) 



v f (t), (3) 



up to corrections of 0(e), where yj = Xf/V6, rj 1 = e5, 
and Vf(t) = v(0,t) is the yet undetermined value of the 
inhibitor at the front position r = 0. 

A dynamical equation for Vf follows from an analysis 
of the outer regions. First we go back to the unstretched 
coordinate system and rescale the spatial coordinate ac- 
cording to y = r/VS. At order unity we find 

vt-VfVy = u + (v)-a 1 v + v vy -a{y + y f )v y + h y < 0, 

(4) 

and a similar equation for y > with u + (v) replaced by 
u—(v). Here u±(v) are the outer solution branches of the 
cubic equation u — u 3 — v = 0. In (4) we assumed the 
field h is constant and took a linear profile J = —ax. For 
ai sufficiently large we may linearize the branches u±(v) 
around v = 0, u±(v) w ±1 — v/2. Inserting (3) in (4) 
and using the approximate forms for u±(v) we find the 
free boundary problem 



Mv = +1 - a(y + Vf)v y + h - 



T]V2 



v(0,t)v y y<0 



Mv = -1 - a(y + y f )v y + h -=v(0,t)v y y > 0, 

?7V2 



v(±oo,t) = v T = TQ 2 H o = [v v ] = 0, (5) 



where M. = d t + q — d£, the square brackets denote 
jumps across the free boundary at y = 0, and v(0, t) and 
y f (t) satisfy (3). 

To solve the free boundary problem (5) we assume the 
system is close to the front bifurcation so that the speed, 
c, of the propagating solutions in the absence of the ex- 
ternal fields is small. We also take the external fields to 
be of order c 3 : a — a n c 3 , h = h n c 3 . We now expand the 
propagating solutions as a power series in c: 



v(y,t,T)=v(°\y) + J2c n v (n) (y,t,T), 



(6) 



n=l 
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where v^°\y), the stationary front solution, is an odd 
function given by v^°'(y) = q~ 2 {e~ qy — 1) for y > 0, and 
T = c 2 t is a slow time scale characterizing the nonsteady 
front motion near the bifurcation. Expanding rj as well, 
V 

in (5) we find 



Vc — c 2 t]i+ c A rj2 + • ■ • , and inserting these expansions 



v. 

where 



(n) 



yy 



1,2,3, (7) 
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The solution of (7) with a zero initial condition (only 
relevant as long as the long fast time asymptotics is con- 
cerned) satisfies the integral equation 
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(9) 



Recall that contains the unknown i/™) evaluated at 
y = 0. Since the origin of the slow time scale is the 
nonsteady front motion we expect «| ™1 to become inde- 
pendent of the fast time scale t as t — ► oo. Substituting 
linit^oo (0, t, T) = (0, T) in (9) and setting y = 
we find a sequence of compatibility conditions. The first, 



for n = 1, is r] c 



2\/2q- 



determines the front bifurcation point. The compatibil- 
ity condition for n = 3 is 



.,(!) 



(0,T) = 



,(1) 



(0,T) 



p) 



(0,T) 3 



(10) 



^7=^- Expressing (10) in terms of v/ 



where 771 = 

v(0,T) w c?;W(0,T) and using c 2 ?yi 



*7 = -^-|(»7c 



4r^ + 3^ + 3 



77 we find 
4 



ft. 



(11) 



where = c 2 vj T . 

Equations (3) and (11) describe the dynamics of fronts 
near a NIB bifurcation, subjected to a constant field h 
and a linearly space dependent J field. In addition to the 



translational degree of freedom, j/j, the dynamics involve 
a second degree of freedom, the value of the inhibitor at 
the front position, Vf, which is responsible for transitions 
between the left and right propagating fronts. Without 
the field J the system (3) and (11) is decoupled and re- 
produces the front bifurcation. The introduction of a 
space dependent field couples the two degrees of freedom 
and affects the front behavior in two significant ways: for 
r\ > rj c (and h ^ 0) it stabilizes a propagating front at 
a fixed position, yf = and for 77 < r\ c it induces 

oscillations between the counterpropagating fronts. The 
frequency of oscillations close to the Hopf bifurcation at 
V = Vc is 
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FIG. 2. (a) Front position, Xf, vs. time for an oscil- 
lating front. The thin line represents the solution to equa- 
tions (3) and (11) and the diamonds are from the numerical 
solution of equations (1). (b) A log- log plot of the oscilla- 
tion frequency, u>, vs. the external field gradient, a. The 
solid line is the relation of equation (12) and the diamonds 
represent numerical solutions of equations (1). Parameters: 
01 = 3.0, e = 0.01, S = 2.77, h = and a = 0.005 in (a). 

To test the validity of equations (3) and (11) we numer- 
ically integrated the original system (1) and compared 
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FIG. 3. The oscillating domain on the left has the same 
frequency as middle domain but the wider separation between 
the J = points produces a larger domain size. The domain 
on the right has the same average width as the middle do- 
main but the larger external field gradient produces higher 
frequency. Parameters: oi = 3.0, e = 0.01, S = 2.77, h = 0. 



oscillating front solutions of (1) with those of (3) and 
(11). The agreement as Fig. 2a shows is very good. In 
Fig. 2b wc plotted the frequency of front oscillations vs. 
the field gradient according to (12) and as obtained from 
(1). Again, the agreement is excellent, and remains good 
even for c of order unity. 

These results suggest various ways to control domain 
patterns. Choosing a periodic J(x) profile, for example, 








20 



40 



60 



80 



FIG. 4. The phase of oscillation is determined by 
the Vf values of the fronts that bound a domain. Op- 
posite signs of Vf give rise to a traveling domain that 
oscillates back and forth (left domain). Equal signs of 
Vf yield a breathing domain (right domain). Parameters: 
oi = 3.0, e = 0.01, S = 2.77, h = 0. 

allows the creation of a periodic pattern of stationary 
(rj > r] c ) or oscillating (i] < r] c ) domains. The period of 
the pattern, the average width of the domains, and the 



frequency and relative phase of oscillations are easily con- 
trollable. Figs. 3 and 4 show the effect of a nonuniform 
triangular profile of J for r] < r/ c . In the absence of J the 
system supports traveling domain patterns. Switching on 
J(x) gives rise to patterns of oscillating domains. Regions 
of J(x) where the gradient a(x) = \ J'{x)\ is steeper yield 
higher oscillation frequencies in accord with (12), while 
wider J triangles yield wider domains (Fig. 3). The rel- 
ative phase of oscillation is controlled by the values of 
Vf for the two fronts that bound a domain. Choosing 
the same sign for Vf gives rise to breathing dynamics 
whereas opposite signs yield back and forth oscillations 
(Fig. 4). For r\ > r\ c arbitrary stationary domain pat- 
terns can be formed with appropriate J{x) profiles; the 
only restriction is the requirement of a minimum domain 
size to guarantee the dominance of J over front inter- 
actions. Similar results are obtained with a nonuniform 
h field and constant J. We expect the main ideas pre- 
sented here to apply to other models exhibiting NIB bi- 
furcations, such as the forced complex Ginzburg-Landau 
equation [4]. 
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